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Abstract 

We study the beam-energy and system-size dependence of (f) meson production (us- 
ing the hadronic decay mode cp K^K~) by comparing the new results from 
Cu+Cu colhsions and previously reported Au+Au collisions at ^snn = 62.4 and 
200 GeV measured in the STAR experiment at RHIC. Data presented are from mid- 
rapidity (|y| <0.5) for 0.4 < px < 5 GeV/c. At a given beam energy, the transverse 
momentum distributions for (p mesons are observed to be similar in yield and shape 
for Cu+Cu and Au+Au colliding systems with similar average numbers of partic- 
ipating nucleons. The (p meson yields in nucleus-nucleus collisions, normalised by 
the average number of participating nucleons, are found to be enhanced relative to 
those from p-\-p collisions with a different trend compared to strange baryons. The 
enhancement for (j) mesons is observed to be higher at -y/sNN = 200 GeV compared 
to 62.4 GeV. These observations for the produced ^(ss) mesons clearly suggest 
that, at these collision energies, the source of enhancement of strange hadrons is 
related to the formation of a dense partonic medium in high energy nucleus-nucleus 
collisions and cannot be alone due to canonical suppression of their production in 
smaller systems. 



Key words: Particle production, Strangeness enhancement. Canonical 
suppression, Quark-Gluon Plasma and Resonances. 
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1 Introduction 



Experimental results from the Relativistic Heavy Ion Collider (RHIC) have 
confirmed the formation of a hot and dense medium in the initial stages of 
high-energy heavy-ion collisions [1]. Thus one of the prerequisites for the for- 
mation of a Quark Gluon Plasma (QGP) [2] in such collisions has been estab- 
lished. High statistics data on 4> meson elliptic flow and yields as a function of 
transverse momentum (px) have been used to support the picture of formation 
of a hot and dense medium with partonic collectivity at RHIC [3] . Evidence of 
(f) mesons being formed by the coalescence of seemingly thermalized ss-quarks 
in central Au-|-Au collisions has also been presented [3]. 

Several interesting features were also observed in the centrality dependence of 
(f) meson production in Au-|-Au collisions at 200 GeV. As one goes from cen- 
tral collisions (average number of participants, (iVpart), > 166) to peripheral 
collisions ((A^part) < 77), the pr spectra showed a gradual evolution from an 
exponential shape to a shape which requires an additional power law type of 
behavior at higher (> 3 GeV/c) [3,4]. At the same time, the average trans- 
verse momentum ((px)) of (f) mesons, dominated by the transverse momentum 
distribution at low px, showed no significant collision centrality dependence 
in Au+Au collisions, unlike what has been seen for other particles of simi- 
lar mass such as anti-protons (p) [4]. The N{(f))/N{K^) ratio was observed 
to be independent of collision centrality in Au-|-Au collisions, in contrast to 
predictions from microscopic transport models like RQMD and UrQMD [5]. 
Both of these results led to the conclusion that (f) meson production may not 
be from KK coalescence and (f) mesons may have decoupled early on in the 
collisions [4]. 

The linear increase of the N{Vt)/N{(f)) ratio with px was proposed as an ob- 
servable to test the recombination picture and hence also provided a test for 
thermalization in heavy- ion collisions [6]. A distinct trend was observed in 
the centrality dependence of this ratio vs. px in Au-|-Au collisions [3]. With 
decreasing centrality, the observed N{fl)/N{(f)) ratio seems to turn over at 
successively lower values of px indicating a smaller contribution from ther- 
mal quark coalescence in more peripheral collisions. Furthermore, in lower 
energy collisions at the SPS [7] and AGS [8], it was observed that the relative 
strangeness production increases with Apart - For similar Apart, the increase was 
found to be slower for larger colliding ions. The possible reason was related to 
variations of space-time density of the participating nucleons and the increase 
in collision density (interactions per fm^) towards the center of the reaction 
volume [7,8]. The measurement of (p production in Cu+Cu collisions, in which 
systems with A^part < 128 are created, is therefore expected to provide more 
precise data to further probe these centrality and colliding ion size dependent 
features. 
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In this letter we report the first results of meson production for rapidities 
\y\ < 0.5 and 0.4 < < 5 GcV/c in Cu+Cu collisions at ^^snn = 62.4 and 
200 GeV. The data were taken by the STAR experiment at RHIC [9]. A 
detailed comparative study of the energy and system size dependence of 
meson production (px spectra, rapidity density and (pt)) is carried out using 
both the Cu+Cu and Au+Au data. 

Several possible mechanisms of (p meson production in nucleus-nucleus colli- 
sions have been reported in the hterature [6,10-12]. Some of these are sup- 
ported by the experimental data [3] which is not true with others [3]. In a 
QGP, thermal s and s quarks can be produced by gluon-gluon interactions [10]. 
These interactions could occur very rapidly and the s-quark abundance would 
equilibriate. During hadronisation, the s and s quarks from the plasma coa- 
lesce to form mesons. Production by this process is not suppressed as per 
the OZI (Okubo-Zweig-Izuka) rule [13]. This, coupled with large abundances 
of strange quarks in the plasma, may lead to a dramatic increase in the pro- 
duction of mesons and other strange hadrons relative to non-QGP p+p col- 
lisions [14] . Alternative ideas of canonical suppression of strangeness in small 
systems as a source of strangeness enhancement in high energy heavy-ion colli- 
sions have been proposed for other strange hadrons (e.g A, S and il) [15]. The 
strangeness conservation laws require the production of an s-quark for each 
s-quark in the strong interaction. The main argument in such canonical mod- 
els is that the energy and space time extensions in smaller systems may not 
be sufficiently large. This leads to a suppression of strange hadron production 
in small coUision systems. These statistical models fit the data reasonably 
well [16]. According to these models, strangeness enhancement in nucleus- 
nucleus collisions, relative to p+p collisions, should increase with the strange 
quark content of the hadrons. This enhancement is predicted to decrease with 
increasing beam energy [17]. So far, discriminating between the two scenarios 
(strange hadron enhancement being due to dense partonic medium formed in 
heavy-ion collisions or due to canonical supression of their production in p+p 
collisions) using the available experimental data has been, to some extent, am- 
biguous. Enhancement of (f){ss) production (zero net strangeness) in Cu-|-Cu 
and Au-|-Au relative to p+p collisions would clearly indicate the formation of a 
dense partonic medium in these collisions. This would then rule out canonical 
suppression effects being the most likely cause for the observed enhancement 
in other strange hadrons [18] in high energy heavy-ion collisions. 



2 Experiment and analysis 

The data presented here were taken at RHIC in 2004 (Au+Au) and 2005 
(Cu+Cu) using the STAR detector [9]. The analysis presented is from the 
data taken by the Time Projection Chamber (TPC) [19]. The TPC magnetic 
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Fig. 1. Upper panel: A typical <j) meson mass peak in Cu+Cu collisions at 200 GeV 
obtained from the K~^K~ invariant mass distribution after subtracting the combi- 
natorial background using mixed events. The distribution is fitted to a Breit-Wigner 
function (solid line) and a linear background function (dashed line) to extract the 
yields. The errors shown are statistical. Lower panel: (j) meson mass peak position 
as a function of pT for various collision centralities in Cu+Cu collisions at 200 GeV. 
Also shown are the results from Monte Carlo calculations for 0-60% centrality using 
embedding techniques (see text for more details) shifted by 50 MeV/c in for 
clarity of presentation. The dashed line corresponds to the PDG value of 1.0194 
GeV/c2 [22]. 
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Table 1 

Collision systems, beam energies, number of events and trigger conditions. 



Collision svstem 


, /,s M AT ( GeV) 


Number of events 




Cu+Cu 


62.4 


8.8 xlO^ 


Minimum Bias 


Cu+Cu 


200 


1.5 xlO'^ 


Minimum Bias 


Au+Au 


62.4 


6.2 xlO^ 


Minimum Bias 


Au+Au 


200 


1.35 xlO^ 


Minimum Bias 


Au+Au 


200 


1.0 xlO^ 


Central Trigger (0-12%) 



field was 0.5 Tesla. Data were taken in both field configurations. The trigger 
conditions and number of events analyzed for different colliding systems at 
^snn = 62.4 GeV and 200 GeV are given in Table 1. The cf) meson spectra for 
Au+Au collisions at 200 GeV using these data sets have been presented else- 
where. [3]. Centrality selection for the Au+Au and Cu+Cu collisions utilized 
the uncorrected charged particle multiplicity for pseudorapidities |?7| < 0.5, 
measured by the TPC. Table 2 shows the (A^part) and (A^bin) values calculated 
using a Glauber model for different centrahties for Cu+Cu coUisions at 62.4 
and 200 GeV and Au+Au collisions at 200 GeV. The corresponding values for 
Au+Au collisions at 62.4 GeV were published previously [20]. 

The (j) meson yield in each pt bin was extracted from the invariant mass (Mj„„) 
distributions of K^K~ candidates after the subtraction of the combinatorial 
background estimated using the event mixing technique [3,4,21]. The charged 
kaons were identified through their ionization energy loss in the TPC. Figure 1 
shows a typical, background subtracted, K^K~ M^nv distribution as obtained 
for 200 GeV Cu+Cu collisions. The resultant distribution is well described by 
a Breit-Wigner function (solid line) plus a linear background function (dashed 

line). The form of the Breit-Wigncr function is = jjj ^^tttftu^ where 

C is the area under the mass peak, V is the full width at half maximum for 
the distribution in GcV/c^ and is the mass of the meson. Figure 1 also 
shows that for pt > 0.7 GeV/c, the mass peak position of the 4> meson agrees 
well with the PDG value of 1.0194 GeV [22]. For < 1-2 GeV/c there is 
a monotonic drop in the value of the fitted mass value with decreasing pt, 
reaching (mass (j) fitted - mass PDG) = -2.5 MeV at pt = 0.5 GeV/c. The 
reconstructed invariant mass distribution of the meson is wider than the 
PDG value (4.26 MeV/c^), decreasing from 9 MeV/c^ to 4.26 MeV/c^ with 
increasing px [23] . The variations in the position of the invariant mass peak 
and its width, at low px, are consistent with the simulation values and are un- 
derstood within the scope of the detector resolution effects [21]. To understand 
these effects, decays to K^K^ and detector response were studied within 
the STAR GEANT framework [24]. The resulting simulated signals were then 
embedded into real events before being processed by the standard STAR event 
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Table 2 

The average numbers of participating nucleons ((A^part)) and binary collisions 
((-^bin)) for various collision centralities in Cu+Cu collisions at ^snn = 62.4 and 
200 GeV and Au+Au collisions at ./s^ = 200 GeV. 



% cs 


200 GeV 


200 G(A' 


200 GeV 


200 G(A' 


(32. 1 G(-V 


62. 1 GeV 


0-10 

10-20 

20-30 

30-40 

40-50 

50-60 

60-70 

70-80 


325.911:1 
234.51?:° 

166.6l^V 
115.51^^ 

76.71^:? 

47.3l^:? 

26.9li:i 

14.111:^ 


939.41^:^ 

590.9li:? 
368.5144:° 
990 1 +35-1 

ZZU.i_34 g 

123.51^^2 
63.9li|^ 
29.51^:1 
12.3118^ 


99.01}:! 

74.6li:o 
53.71^:° 

37.81°:,^ 
26.2l°:l 

17.21°:^ 


188.8111:^ 
123.6l|^ 
77.611:^ 
47.7l|? 

29.2ii:^ 
16.81°:? 


96.41^:^ 
72.21°:^ 

51.81°:! 

36.21°:^ 
24.91°:^ 

16.31°:^ 


161.81^^5^ 
107.5l°:i 

68.4l|^ 

25.9li:° 

15.11^:^ 



reconstruction. These data were then processed hke real data and analyzed to 
reconstruct the embedded [3,4,21,23]. Embedding simulations were also used 
to obtain the (p meson acceptance and reconstruction efficiency [21,23]. The 
product of the acceptance and reconstruction efficiency was found to increase 
from 3% at px = 0.5 GeV/c to about 40% at = 3 GeV/c for central Cu+Cu 
collisions. The centrality dependence of these values were found to be small 
for Cu+Cu collisions. At higher (3-5 GeV/c), the efficiency was found to 
remain constant. The other important corrections applied to the data were re- 
lated to the vertex finding efficiency which was ~ 92.5% and the correction for 
branching ratio of 49.2% for the channel cj) K'^K". A more detailed descrip- 
tion of the meson mass peak position, width of the meson invariant mass 
distribution, variation of the reconstruction efficiency with collision centrality 
and Pt, and the general procedure for obtaining the signal and constructing 
mixed events are discussed elsewhere [23]. 

Systematic errors for the meson spectral measurements in Cu+Cu collisions 
include uncertainties from the following sources: Uncertainties in (p meson 
reconstruction efficiency (~ 8-14%), Kaon identification from dE/dx (8%), 
Kaon energy loss corrections (~ 3-4%), Residual background shape (4%) and 
magnetic field configuration (~ 3%). The systematic errors from all the above 
sources have been added in quadrature. Systematic errors for the meson 
spectra are similar at both energies (62.4 and 200 GeV). The total systematic 
errors for yields at both energies are estimated to be ^ 18% over the entire px 
range studied. A discussion on systematic errors for Au+Au collisions, dN/dy, 
and (px) can be found in Ref. [3,4,23]. 
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Fig. 2. Midrapidity (|y| < 0.5) transverse momentum spectra of (/> mesons for various 
collision centrality classes for Cu+Cu collisions at ^snn = 62.4 and 200 GeV. 
To study the system size dependence, comparison of 40-50% Au+Au spectra to 
10-20% Cu+Cu spectra at 200 GeV, and 40-60% Au+Au spectra to 20-30% Cu+Cu 
spectra at 62.4 GeV are shown. These centralities for the two colliding systems have 
similar (A'part) values as outlined in Table 2. The errors represent the statistical and 
systematic errors added in quadrature. They are found to be within the symbol size. 
The spectra are fitted to a Levy function discussed in the text. 

3 Transverse momentum distributions and yields 



Figure 2 shows the (j) meson yields from Cu+Cu collisions at 62.4 and 200 GeV 
for 0.4 < pt < 5 GeV/c and various collision centralities. The spectra are well 
described by a Levy function of the form ^ — '^^^ , = j—j , „ — r;-, where 

rriT = \Jp\ + m?^. A, Ti^evy, and n are the parameters of the function. In the 
limiting case of 1/n — >^ 0, the Levy distribution approaches an exponential 
function. The parameters TLevy and n have similar values for the Cu+Cu and 
Au+Au systems with similar (A/'part) at 200 GeV. This reflects the similar 
shape for the meson spectra in both collision systems at a given energy and 
(A^part)- A comparison of mesons spectra for 40-50% central Au+Au ((Apart) 
= 76.7) and 10-20% central Cu+Cu ((A^part) = 74.6) collisions at 200 GeV is 
shown in Fig. 2 (left panel). Similar results for 40-60% central Au+Au ((Apart) 
= 59.9) and 20-30% central Cu+Cu ((A^part) = 51.8) collisions at 62.4 GeV are 
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Table 3 

dN/dy and (px) for ^ mesons produced in Cu+Cu collisions at ^/sNN = 200 and 
62.4 GeV for various collision centralities. The errors include both systematic and 
statistical errors added in quadrature. 



% centrality 


dN/dy 


dN/dy 


< PT > (MeV/c) 


< PT > (MeV/c) 




200 GeV 


62.4 GeV 


200 GeV 


62.4 GeV 


0-10 


2.3 ± 0.3 


1.3 ± 0.2 


935 ± 60 


881 ± 61 


10-20 


1.6 ± 0.2 


0.8 ± 0.1 


901 ± 64 


857 ± 65 


20-30 


1.1 ± 0.2 


0.6 ± 0.1 


897 ± 62 


848 ± 57 


30-40 


0.7 ± 0.1 


0.4 ± 0.1 


885 ± 60 


835 ± 57 


40-50 


0.4 ± 0.1 


0.24 ± 0.04 


869 ± 59 


793 ± 57 


50-60 


0.26 ± 0.05 


0.14 ± 0.03 


852 ± 58 


771 ± 56 



also shown in the same figure on the right panel. The ratios of the (f) meson 
Pt spectra for Au-|-Au and Cu-|-Cu systems with similar (A^part) agree within 
~ 10%. This is further quantified by studying their rapidity density {dN/dy) 
and (pt) for both colliding systems. 

Figure 3 shows dN/dy, dN / dy / {N^g^j.^) and < > as a function of (iVpart) 
for Cu-|-Cu and Au-|-Au collisions at 62.4 and 200 GeV. Results from p-\-p 
at 200 GeV and 63 GeV, obtained from the STAR [4] and ISR [25] exper- 
iments respectively, arc also included for comparison. At 63 GcV the da/dy 
for mesons at < y < 0.33 was reported to be 0.44 ±0.11 (sys) ±0.1 
(stat) mb. These data, together with values of 36 and 42 mb for p+p inelastic 
cross-sections at 63 and 200 GeV respectively, have been used to get the cor- 
responding dN/ dy values shown in the figure. The dN/dy and < px > values 
as obtained for the Cu±Cu collisions are also presented in Table 3. Both at 
62.4 and 200 GeV, all three quantities viz dN/dy, dN/dy / (Nps^j-t) and (px) 
scale with (Apart)- These findings seem to indicate that the general features 
of (f) meson production characterized in terms of dN/dy and (px) at a given 
energy (62.4 or 200 GeV) do not depend on the colliding ion species studied, 
but depend on the (A'part) of the collision. It will be interesting to see whether 
the same is true for other produced hadrons at RHIC. However, for a given 
(Apart), both dN/dy and (px) are observed to be lower for 62.4 GeV when 
compared to 200 GeV. This is in contrast to what has been seen at lower 
energies at AGS and SPS with smaller colliding systems [7,8]. At those lower 
energies, for similar (Apart), the strange hadron production was higher while 
at RHIC, due to higher center of mass energy, a hotter and denser medium 
is expected to form with a very low net baryon density at midrapidity [1], 
leading to the observed differences. 
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Fig. 3. Upper panels: dN/dy at midrapidity for (j) mesons for various collision cen- 
trality classes in Cu+Cu and Au+Au at ^snn = 200 GeV and 62.4 GeV. Also 
shown are the results from p+p collisions. Middle panels: same as above, but for 
dN / dy / {Nps,T-t) . Lower panels: Average transverse momentum {{pt)) for mesons 
at midrapidity for various event centrality classes for Cu+Cu and Au+Au collisions 
at y^SNN = 62.4 GeV and 200 GeV. The (px) for 4> niesons in p+p collisions are 
also shown. The error bars represent the statistical and systematic errors added in 
quadrature. 

4 Nuclear modification factor 



Now we look at the dependences of the nuclear modification factor, for the 
(j) meson, both in terms of A^part and A^bin. For A^part, this factor is given by 

Popart/ N ^ d^NAA/PTdydpT/ (A/"part) 
^AA [PT) d^a^Jdydpr/a^l^^ ' 

To get the corresponding i?^^'"(pT), one needs to replace (A^part) by (A^bin) 
in the above expression. The results, as shown in Fig. 2 and Fig. 3 would 
lead to very similar results on -Raa ' both Cu+Cu and Au+Au systems 
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Fig. 4. Upper panels: A'^part scaled (-Raa'") nuclear modification factor as a function 
of of mesons for 0-10% and 20-30% Cu+Cu and Au+Au collisions at ^snn 
= 200 GeV. Lower panel: Same as above for A^bin scaled (-Raa") iiuclear modifica- 
tion factor. The error bars represent the statistical and systematic errors added in 
quadrature. The shaded band in upper panel around 1 at pT = 4.5 — 5.5 GeV/c 
in the right side reflects the uncertainty in (A'^part) and that on the lower panel 
for (A^bin) calculation for central Au+Au collisions. The respective uncertainties for 
central Cu+Cu collisions are of similar order. 

for collisions having similar (A^part)- In view of this, we only present here a 
comparison of the nuclear modification factors (in terms of Nun and Npart) of 
i?^A" Raa^^ for Cu+Cu and Au+Au collisions. For such a comparison only 
centralities corresponding to similar fraction of total hadronic cross-section 
were considered. The Raa for 4> mesons in 200 GeV Cu+Cu and Au+Au 
collisions for 0-10% and 20-30% collision centralities (up to = 4 GeV/c) at 
200 GeV are shown in Fig. 4. 

Within the errors, the -Raa"^' values for 0-10% central Cu+Cu and Au+Au 
collisions at 200 GeV are seen to be similar in shape and yields. However, 
for 20-30% collisions and at other collision centralities (which are not shown 
in the figure) the Au+Au results are higher than Cu+Cu results for most of 
the pt range studied. The results for the central most Cu+Cu and Au+Au 
collisions studied are consistent with the observation that dN/dy/ (A^part) and 
< Pt > are constant as a function of (A/'part) for (A^part) > 90 (Fig. 3). 

At the same collision centralities, the ratio (-^biii^^)/(^biii*"") about ~ 1.5 
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times larger than the ratio (N^^^"^) / {N^^^"") . This is reflected in the T^aa"- 
As one can see from Fig. 4, -Raa" 0-10% Cu+Cu is higher than that of 
Au+Au coUisions, for < 3 GeV/c. Both the modification factors at > 3.5 
GeV/ c are below unity, showing the characteristics of parton energy loss in hot 
and dense medium formed in central heavy-ion collisions. For 20-30% central 
collisions, the similarity between i?AA° for Cu-|-Cu and Au-|-Au collisions seems 
to extend to lower px (~ 1-5 GeV/c). It may be interesting to use the nuclear 
modification factor of mesons to investigate the differences in energy loss of 
quarks and gluons in the medium formed in heavy- ion collisions [26]. This is 
because (f) mesons in central collisions are formed from coalescence of s and s 
quarks [3], which presumably are formed by gluon interactions in the initial 
stages of the coUision. 



5 meson production and strangeness enhancement 



The ratio of strange hadron production normalized to (A^part) in nucleus- 
nucleus collisions relative to corresponding results from p+p collisions at 200 
GeV are shown in the upper panel of Fig. 5. The results are plotted as a func- 
tion of (A^part)- K~ [27], A and S + S [18] are seen to show an enhancement 
(value > 1) that increases with the number of strange valence quarks. Fur- 
thermore, the observed enhancement in these open-strange hadrons increases 
with collision centrality, reaching a maximum for the most central collisions. 
However, the enhancement of (f> meson production from Cu-|-Cu and Au-|-Au 
collisions shows a deviation in ordering in terms of the number of strange con- 
stituent quarks. More explicitly, this enhancement is larger than for K~ and 
A, at the same time being smaller than in case of S -|- S. Despite being differ- 
ent particle types (meson-baryon) and having different masses, the results for 
K~ and A are very similar in the entire centrality region studied. This rules 
out a baryon-meson effect as being the reason for the deviation of mesons 
from the number of strange quark ordering seen in Fig. 5 (upper panel). The 
observed deviation is also not a mass effect as the enhancement in meson 
production is larger than that in A (which has mass close to that of the 0) . 

In heavy-ion collisions, the production of mesons is not canonically sup- 
pressed due to its ss structure. In low energy p+p collisions at ^/s = 3.6 
GeV, meson production is suppressed due to the OZI rule [28]. In p+p col- 
lisions at ^/s — 6.84 GeV violation of this rule has been reported [29]. At this 
higher energies production through channels accompanied by non-strange 
hadrons was found to dominate strongly over its production in channels ac- 
companied with strange hadrons. Measurements of production in proton- 
nucleus collisions at ^/sNN = 27.4 GeV have also shown that it takes place 
primarily by other than OZI allowed processes [30]. Experiments studying 
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Fig. 5. Upper panel: The ratio of the yields of K~ , (p, A and H + H normalized to 
impart) in nucleus- nucleus collisions to corresponding yields in inelastic p+p colli- 
sions as a function of (iVpart) at 200 GeV. Lower panel: Same as above for (p mesons 
in Cu-|-Cu collisions at 200 and 62.4 GeV. The p+p collision data at 200 GeV are 
from STAR [4] and at 62.4 GeV from ISR [25]. The error bars shown represent the 
statistical and systematic errors added in quadrature. 



inclusive production off protons by hadrons at incident momenta 63 and 
93 GeV/c also show that the production of mesons are from OZI allowed 
processes [31]. Experiments on the production of mesons in p+p collisions 
near threshold have shown a large enhancement of the cross section ratio 
a{pp — >• pp(l))/a{pp — i> ppu) [32] compared to that predicted by the OZI 
rule [33]. This ratio is sensitive to the basic feature of the rule, which states 
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that proceses with disconnected quark hnes between initial and final states are 
suppressed compared to those where the incident quarks continue through to 
the exit channel. The p+p collisions at RHIC are at an energy which is ~ 25 
times higher than energies where violations of the OZI rule were reported [29] . 
The meson enhancement in heavy-ion collisions shows an increasing trend 
with centrality (Fig. 5). Prom this, we conclude that the observed enhancement 
of (p production in heavy-ion collisions may not be due to OZI suppression of 
4> production in p+p collisions. 

The observed enhancement of (j) meson production then is a clear indication for 
the formation of a dense partonic medium being responsible for the strangeness 
enhancement in Au-|-Au collisions at 200 GcV. Furthermore, (p mesons do not 
follow the strange quark ordering as expected in the canonical picture for the 
production of other strange hadrons. The observed enhancement in meson 
production being related to medium density is further supported by the energy 
dependence shown in the lower panel of Fig. 5 . The (p meson production rela- 
tive to p+p collisions is larger at higher beam energy, a trend opposite to that 
predicted in canonical models for other strange hadrons. Earlier measurements 
have indicated that (f) meson production is not from coalescence of KK and 
minimally affected by re-scattering effects in the medium [4] . Recent measure- 
ments indicate that (p mesons are formed from the coalescence of seemingly 
thermalized strange quarks [3]. All these observations put together along with 
the observed centrality and energy dependence of (p meson production (shown 
in Fig. 5) indicate the formation of a dense partonic medium in heavy-ion 
collisions where strange quark production is enhanced (possible mechanisms 
could be as discussed in Refs. [10,14]). This in turn suggests that the observed 
centrality dependence of the enhancement for other strange hadrons (shown in 
Fig. 5) is likely to be related to the same reasons as in the case of the (p meson, 
that it is due to the formation of a dense partonic medium in the collisions. 
These experimental data rule out the possibility of canonical suppression be- 
ing the only source of the observed strangeness enhancement at beam energies 
of 200 GeV. 



6 Summciry 

We have presented a study of the energy and system size dependence of (p 
meson production using the p+p, Cu+Cu and Au+Au data at a/snn = 62.4 
and 200 GeV. The px spectra are measured at midrapidity {\y\ < 0.5) over the 
range 0.4 < p^ < 5 GeV/c. These measurements provide new experimental 
results showing that at a given beam energy the transverse momentum spectra 
in both shape ((pt)) and yields (dN/dy) are similar in Cu+Cu and Au+Au 
for collisions with similar (iVpart)- In addition to observing similarity in the 
(p meson distributions for Cu+Cu and Au+Au collisions with similar (A'part), 
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the (iVpart) scaled nuclear modification factors are observed to be similar for 
the 0-10% central Cu+Cu and Au+Au collisions at 200 GeV. However, such a 
similarity is not seen for other collision centralities. The corresponding results 
for the nuclear modification factor, scaled by the number of binary collisions, 
are in general found to be higher for Cu+Cu compared to Au+Au collisions. 



The enhancement in the meson production has been studied through the 
ratio of the yields normalized to (Npg^-i) in nucleus-nucleus collisions to corre- 
sponding yields in p+p collisions as a function of (A^part)- The centrality and 
energy dependence of the enhancement in meson production clearly reflects 
the enhanced production of s-quarks in a dense medium formed in high en- 
ergy heavy-ion collisions. This then indicates that the observed enhancements 
in other strange hadron {K^ .A and S + S) production in the same collision 
system are likely to be due to the similar effects and not only due to canonical 
suppression of strangeness production. At RHIC the colliding beam energy 
is high, so it is very unlikely that the observed enhancement in heavy-ion 
collisions is due to OZI suppression of (p production in p+p collisions. 



The enhancement in the meson production deviates from the number of va- 
lence s-quark dependence observed for other strange hadrons. The results from 
mesons lie in between those from single valence s-quark carrying hadrons K~ 
and A, and double valence s-quark carrying hadrons S + S. Comparisons with 
other strange hadrons rule out the possibility of this being a baryon-meson or 
mass effect. The exact reason for the observed deviation of the enhancement 
factor for the meson from the valence strange quark dependence observed for 
other strange hadrons is not clear. It could be due to the effect of light-flavor 
valence quarks in the other strange hadrons or due to the net strangeness 
being zero in (p mesons. 
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